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1.0 INTRODUCTION 


Previous reports reviewed the earlier efforts to disperse fog [1] 
and examined the theoretical aspects of the charged particle fog dis- 
persion technique in great detail [2]. The charged particle concept 
appeared to be a viable technique based on the review of this previous 
work, but it was concluded that further experimental investigation would 
be required to remove remaining doubts. This report details the design 
of a prototype charged particle fog dispersal unit and outlines the 
experimental program to evaluate the unit. The construction and evalua- 
tion of the unit will be performed in the next phase of this work. The 
techniques used to design the unit are discussed in Section 2.0, and 
Section 3.0 contains detailed plans for the prototype unit. The appen- 
dix contains a bibliography of work performed in electrofluid dynamic 
power generation that was found to be helpful in the design of the 
prototype unit. 


2.0 DESIGN OF PROTOTYPE CHARGED PARTICLE FOG DISPERSAL UNIT 


2.1 Previous Designs 

Previously used charged particle fog dispersal units and their 
operation were described in the earlier reports [1,2] and in References 
[3] through [10]. This previous work will not be reviewed again here, 
but the results will be freely used. During the course of the design of 
the prototype unit the experience that had been gained by other organi- 
zations, i.e., the Curtiss-Wright Corporation, Marks Polarized Corpora- 
tion, University of Dayton Research Institute, and the Aerospace 
Research Laboratories at Wright-Patterson Air Force Base, on the design 
of electrofluid dynamic (EFD) energy conversion devices was found to be 
of considerable assistance because of the great similarity of EFD devices 
and charged particle fog dispersion devices. A bibliography of work 
related to EFD devices is included in the appendix. Certain features of 
the nozzles were found to have significant influences upon increasing 
the amount of power generated by EFD devices. These features have been 
incorporated into the design of the prototype fog dispersal unit and 
flexibility has been included to determine whether improved performance 
can be obtained by utilizing these features. 

2.2 Description of Prototype -Unit 

The features of the prototype charged particle fog dispersal unit 
to be described in this report are shown in Table 1. Previous investi- 
gations of fog dispersal systems [2] or electrofluid dynamic energy 
conversion systems (see appendix) have used various types of charge 
carriers (ions [3,11]; water droplets [4,5,6]; and helium bubbles [7]) 
that were charged either by electrostatic induction [8] or by charge 
transfer in a corona discharge [12,13,14,15]. The charge carriers were 
accelerated by electrostatic means [16] or by viscous coupling to an 
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TABLE 1. Features of Prototype Nozzle. 


• Moist air accelerated in supersonic nozzle having fixed exit Mach 
number. 

• Water droplets formed by homogeneous nucleation in expanding nozzle 
flow. 

• Water droplets charged in corona discharge inside nozzle. 

• Variable nozzle length to sonic radius, 
t Variable corona discharge location. 

• Adjustable corona needle location. 

• Variable stagnation chamber specific humidity. 

• Variable corona needle shape. 

• Constructed from a transparent material. 


accelerating airstream [17]. It was concluded from this previous work 
that the optimum efficiency of charge transport could be obtained by 
charging small water droplets in a corona discharge and accelerating 
these by viscous coupling with air in a supersonic nozzle [12,17]. 

The prototype charged particle fog dispersal unit operation will 
now be described in detail. Moist air is expanded through the sonic 
throat and expansion nozzle. At some point, condensation from sponta- 
neous (homogeneous) nucleation occurs. Heterogeneous condensation is 
not usually considered an effective mechanism for nozzle flows [18]. In 
moist air expansions the nucleation and growth processes are uncoupled 
since homogeneous nucleation is over before any appreciable amount of 
growth has occurred [19]. The droplets grow quickly for some distance 
through the "condensation zone" and then more slowly up to the nozzle 
exit [19,20]. The location of the point where homogeneous nucleation 
takes place is dependent upon the average temperature gradient along the 
nozzle because, for a fixed exit Mach number, the amount of supercooling 
increases with an increase in the average temperature gradient along the 
nozzle (dT/dx) . This temperature gradient will be varied by varying the 
length of the nozzle. 
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The location for homogeneous nucleation and the droplet growth rate 
are extremely sensitive to the relative humidity in the stagnation region 
and, in addition, the droplet growth rate in the "condensation zone" 
depends on the length of the nozzle [18]. In general, small drops are 
produced in short nozzles and large drops in long nozzles. At high 
specific humidities it has been noted that the entire flow is foggy down- 
stream of the condensation location. Lawson [21] found that the maximum 
current for electrofluid dynamic energy conversion nozzles having an 
exit Mach number of 2.0 was obtained with a stagnation chamber relative 
humidity of 17 percent. This rather low value of relative humidity 
required to obtain optimum results will be explained in Section 2.3. 

The corona discharge must be placed after condensation has occurred 
and also at the location where the droplets have grown to a desirable 
size to be able to carry the necessary charge [2]. Therefore, several 
corona locations will be examined. In addition, several corona shapes 
will be examined and the needle location will be adjustable relative to 
each corona discharge location. 

The nozzle will be made out of a transparent plastic in order to 
visualize the flow and corona discharge within the nozzle. Further 
explanation of the features of the prototype nozzle will be given later. 

A schematic of the prototype fog dispersal nozzle is shown in 
Figure 1. Atmospheric air will be compressed to the proper stagnation 
pressure and then filtered to remove all moisture, oil droplets, and 
solid particles. The air then will be mixed with a controlled amount of 
atomized water--to control the specific humidity--in a mixing tube. The 
purpose of the mixing tube is to allow for complete evaporation of the 
water into the dry air. The swirl induced by the mixing tube will be 
removed by a honeycomb. The humid air will be decelerated into the 
stagnation chamber and then expanded through the converging/diverging 
nozzle until it will emerge at the nozzle exit at a supersonic Mach 
number. The diverging portion of the nozzle will be conical with the 
length much greater than the diameter at the throat. Therefore, the 
flow in the nozzle will be very nearly cylindrical. 
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Figure 1 Prototype fog dispersal nozzle. 


The water vapor will condense near the throat during the expansion 
process. The droplets that will be formed will be small enough to be 
carried by viscous drag along with the air. Thus, they will be continu- 
ously accelerated in the nozzle. 

The corona needle will be attached to a veined support in the 
stilling chamber. The attractors can be placed at two locations and two 
lengths of nozzle will be examined. The overall length of the nozzle 
portion of the unit will be approximately 31.75 cm (12.5 in) to 40.64 cm 
(16 in) and 7.62 cm (3 in) in diameter. The entire system can be disas- 
sembled quickly and reassembled in a new configuration. 

The estimated characteristics of the nozzle are shown in Table 2. 

The exit Mach number is 1.35. This is the Mach number that has been 

used by Gourdine [5] and the choice of the Mach number will be explained 

-3 3 

later. A 7.08 x 10 m /s (15 SCFM) compressor will be used which for 
an exit Mach number of 1.35 fixes the maximum diameter of the nozzle 
throat at 0.43 cm (0.168 in). The maximum water flow rate into air that 
is initially dry will be 1.42 x 10"^ m^/s (0.135 gal/hr), which will 


result in a relative humidity in the 
The stagnation pressure will be 3.01 

TABLE 2. Nozzle 

stagnation region of 100 percent. 
X 10^ N/m^ (43.6 psia) for an 

Characteristics. 

Exit Mach number 

1.35 

Minimum diameter 

0.43 cm (0.168 in) 

Airflow rate 

7.08 X 10"^ m^/s (15 SCFM) 

Maximum water flow rate 

1.42 X 10'^ m^/s (0.135 gal/hr) 

Stagnation pressure 

3.01 X 10^ N/m^ (43.6 psia) 

Nozzle length to sonic 


diameter ratio 

40/60 

Mach number at corona 

1.07/1.15 

Current (maximum) 

52 X 10“ amps 

Charge to mass ratio 

0.37 coul/kg 
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atmospheric pressure of 1.01 x 10 N/m (14.7 psia). The stagnation 
pressure is fixed by a desire to eliminate shocks at the exit of the 
nozzle since shocks will result in breakup and evaporation of the water 
droplets. The nozzle length to diameter ratio will be 40 or 60 based on 
the diameter at the nozzle throat. The corona region will be placed 
either at a Mach number of 1.07 or 1.15. The maximum estimated current 
will be 52 microamps and that will result in a charge to mass ratio of 
0.37 coul/kg. The exit Mach number will later be shown to be a com- 
promise between a desire to increase the average temperature gradient 
along with nozzle, to increase the amount of supersaturation, and the 
need to keep the exit Mach number low in order to have as high a charge 
density as possible without a breakdown field at the exit of the nozzle. 

The electrical connections to the nozzle are shown in Figure 2. 

The corona needle will be grounded whereas the attractor will be placed 
at a high but variable positive potential. This will result in a 
negative corona which is the desirable type of corona because of the 
presence of electronegative gases in air [2,22,23,24]. Microammeters 
will be used to measure the cathode or needle current and the attractor 
or anode current. The properties of the discharge will be adjusted 
until the attractor current is nearly zero and that will result in the 
maximum jet current. 

The nozzle controls are shown schematically in Figure 3. The air 
from the compressor will be filtered to remove oil and water droplets. 

It will then be divided; part of the air being used to pressurize the 
water storage bottle and part of the air going directly into the nozzle. 
Water leaving the storage bottle will be strained and the flow rate will 
be measured with a rotameter. The water will flow through a capillary 
tube that will atomize the water when it meets with the air. The 
air/water mixture will then go through the mixing tube where the water 
will evaporate before reaching the nozzle. All quantities, i.e., air 
pressure, water pressure, and water flow rate, can be continuously 
controlled in order to optimize the performance of the nozzle. 
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Figure 3 Nozzle controls. 



2.3 Design of Prototype Unit 


In this section the rationale for the design decisions that were 
described in the previous section will be given. First, the variable 
needle location is examined. Studies on electrofluid dynamic power 
generation devices indicated that the output current can be increased by 
adjusting the position of the needle relative to the attractor [25]. 

The apparatus used in these studies is shown in Figure 4. High-pressure 
air that proceeded through the corona discharge located at the nozzle 
sonic orifice was mixed with medium-pressure secondary air before pro- 
ceeding through a grounded electrode and into the conversion duct. The 
configuration is shown to indicate the relative position of the corona 
needle to the attractor nozzle. By varying the needle location relative 
to the attractor nozzle, the current is found to have a maximum--as 
illustrated in Figure 5. Several corona needle locations relative to 
the attractor will be examined during testing of the nozzle. 

It should also be noticed that the results shown in Figure 5 indi- 
cate that improved performance was obtained with a nail-head needle (a 
#8 nail was also used [21]) compared to a sharp needle. It was felt 
that the low-pressure region and wake that were created downstream of 
the nail head was advantageous for charging the droplets because the 
droplets were trapped for a period of time in the separated region where 
the corona discharge current is increased and therefore improved perfor- 
mance was obtained. Figure 6 indicates the geometry of the corona 
needles that will be examined in the study to try to obtain the maximum 
current from the nozzle. 

It is desirable to estimate the characteristics of the droplets 
that can be expected to be leaving the nozzle. These estimated charac- 
teristics are shown in Table 3. The following assumptions were used in 
making these estimates. The water flow rate was assumed to be 1.42 x 
10"^ m^/s (0.135 gal/hr) which would give a relative humidity of 100 
percent in the stagnation chamber at a temperature of 21.1° C (70° F). 

It was assumed that all of the droplets were charged to the breakdown 
field at their surface. This then allowed the number of electrical 
charges per particle and the mobility of the droplet to be estimated [2]. 
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Figure 4 Typical electrode configuration of EFD test rig [25]. 
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Figure 6 Corona needles. 


TABLE 3. Characteristics of Nozzle Droplets. 


Water flow rate = 1.42 x 10"^ n?/s (0.135 gal/hr) 

All particles charged to breakdown field at surface. 


Particle 
Radi us 
(m) 

Particle 

Flux 

(#/sec) 

Maximum 

Charges 

Per 

Particle 

(#) 

Maximum 

Current 

(amp) 

Mobi 1 i ty 
(m2/volt sec) 

Charge 

Density 

(coul/m3) 

5x10"^ 

2.7x10^^ 

0.05 

2.25x10"^ 

1.10x10'^ 

2.43x10"^ 

10'^ 

3.4x10^^ 

0.21 

1.1 2x1 0‘^ 

1.12x10"^ 

1.21x10'^ 

5x10"^ 

2.7x10^^ 

5.20 

2.25x10"^ 

1.41x10"^ 

2.43x10"^ 

10"^ 

3.4x10^^ 

20.80 

1.12x10'^ 

1.82x10"^ 

1.21x10'^ 


Because of the Cunningham correction factor [2], the mobility is seen to 

be approximately constant for particles in this size range and is approxi- 
-7 2 

mately 10 m /volt sec. 

If the flow within the nozzle is in equilibrium, then the equilib- 
rium drop radius can be estimated from a knowledge of the relative 
humidity at any particular point along the nozzle. This estimate indi- 
cates that the droplets within the nozzle are of extremely small diameter 

[18,19,20]. The droplets usually have diameters of 10 A° to 400 A° (1 A° 
-4 \ 

= 10 microns). Therefore, the radii of the particles examined in 
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Table 3 are quite small, and it is felt that this is the probable range 
for the droplet size from these nozzles [5,19]. Thus, knowing the water 
flow rate, the particle flux rate can be determined, assuming spherical 
particles. From the maximum charge per particle and the particle flux 
rate, the maximum current from the nozzle can be determined. 

While the calculations that were performed to obtain the results in 
Table 3 are straightforward, the results are in all probability incorrect 
because no use has been made of the equations of electrodynamics. 

Magnetic field effects are negligible because of the small currents and 
only electrostatic effects need to be considered. Therefore, the elec- 
tric field is described by the Poisson equation. Minardi [26] computed 
the performance characteristics of electrofluid dynamic power generation 
nozzles assuming that the charge carrier was a colloid such as a water 
droplet. He assumed that he had a very long cylinder through which the 
gas and charged colloids flowed, that the gas velocity was constant 
within the long cylinder, that the charge density was constant within 
the cylinder, and that there was no slip velocity between the colloid 
and the airstream. The last assumption amounts to assuming that the 
mobility of the colloids is zero. Turbulent spreading and diffusion of 
the charged particles were neglected. The Poisson equation was then 
used to determine the electric field within and external to the nozzle 
[27,28,29]. The maximum field was found to occur at the nozzle exit. 

The parameters used in his analysis are the following: 

PghG 

D = — F — = normalized charge density 

^o^bo 


G = 




= configuration factor if R/h << 1 


0.833 + 0.167M? 

where is the electric charge density, R is the radius of a channel of 
length h, is the permittivity of the gas, is the gas breakdown 
electric field in the stagnation chamber, M is the gas Mach number, and 
M* is the Mach number based on the local speed of sound. 
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The factor G is shown in Figure 7. For a very long channel, where 
the radius is much less than the length of the channel, the factor G can 
be approximated by the equation given and is then a very small number. 
Lines of constant normalized output power from the electrofluid dynamic 
nozzle are shown in Figure 8 as a function of M* and V", the normalized 
potential applied to the nozzle. The results of this figure can be used 
in the prototype nozzle because the line V = 0 gives the values of D for 
breakdown in the channel as a result of the charge density alone. 


G 


Ax i symmetric 



R/h or w/h 

Figure 7 Configuration factor of G for both the two-dimensional and 
axisymmetric geometries [26] . 


As previously mentioned, the maximum electric field occurs at the 
nozzle exit and therefore breakdown will occur there first. For a given 
nozzle Mach number, M* can be determined. Then from Figure 8 the value 
of D can be obtained. It should be noted that lines of constant D 
intersect the M* axis at V = 0. From the value of D obtained from this 
figure, the maximum value of the charge density within the nozzle can 
be determined and then the maximum current is given by 

I . = p VttR2 

3 e 

where V is the gas velocity. 
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Figure 8 Lines of constant normalized power as a function of M* and V 
for a gas with k = 1.4. The field strength at the generator 
entrance is [26]. 


The above analysis was applied to the prototype nozzle assuming 
that the nozzle was a cylinder with radius equal to the radius at the 
nozzle throat. The breakdown field was found to occur at the nozzle 
exit when the charge density reached 5.64 x 10" coul/m (M* = 1.266 for 
an exit Mach number of 1.35) for both nozzle length to diameter ratios 
that will be used in the prototype nozzle. This yields a maximum proto- 
type nozzle current of 52.3 microamps, which is considerably less than 
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the estimated currents shown in Table 3. This maximum nozzle current 

would occur for any particle radius within the range considered and it 

therefore appears that the nozzle current is controlled by the breakdown 

of the field at the nozzle exit rather than the details of the charging 

process, provided it is sufficiently efficient. In addition, under the 

assumptions that were made in calculating Table 3, it is seen that the 

stagnation chamber relative humidity can be greatly reduced and the 

_g 

nozzle will deliver the same current (2.3 percent for 5 x 10 m radius 
droplets to 47 percent for 10"^ m radius droplets), i.e., the water flow 
rate can be greatly reduced. Lawson obtained optimum results with a 
relative humidity of 17 percent [21]. 

As a verification of the applicability of this analysis to the 
current from a fog dispersal nozzle, it was applied to the nozzle used 
by Gourdine [5]. The analysis estimates that the maximum current would 
be 99 microamps while the current was measured to be 100 to 135 microamps. 
Therefore, it is felt that the analysis of Minardi can be applied to the 
prototype nozzle and will give a reasonably accurate estimate of the 
maximum current that can be expected from the nozzle. During testing, 
the various features of the prototype nozzle will be studied to determine 
the maximum current the nozzle will deliver. 

Notice that Figure 8 indicates that the value of D for breakdown of 
the field at the exit decreases as M* increases or as the Mach number of 
the nozzle exit increases. Therefore, increased charge density can be 
obtained by lowering the Mach number. But the condensation process 
requires a minimum value of the mean temperature gradient across the 
nozzle, which decreases as the exit Mach number is increased (for a 
fixed length to sonic diameter ratio). Gourdine [5] used nozzles having 
an exit Mach number of 1.35 and it is proposed that the prototype nozzle 
that is described in this report have the same Mach number since this 
appears to be a good compromise between these two competing effects. 

It was desired to choose a size for the prototype nozzle that could 
be economically manufactured in large numbers and would also place the 
required number of charges into the fog in a reasonable length of time. 

The design trade-offs that are possible can be obtained by examining the 
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scaling laws for the corona discharge [12,22]. These scaling laws, 
which are summarized in Table 4, indicate the effect of varying the 
nozzle dimensions and the stagnation pressure upon the discharge. If 
the dimensions of the nozzles are doubled, then the current can be 
increased by a factor of two. However, the nozzle area then has been 
increased by a factor of four which results in a four-fold increase in 
the air volume flow rate in standard cubic feet per minute. This then 
requires a compressor of four times capacity. However, that compressor 
is approximately four times as expensive as the one for the original 
nozzle while leading to only a doubling of the nozzle current. There- 
fore, nozzles of a moderate size are more economical than half the 
number of nozzles that are twice as large since the compressor is the 
most expensive component. Gourdine suggested that the medium-sized 

nozzles that he had examined would be suitable for an airport application 

-3 3 

[30,31]. Because of the availability of a reliable 7.08 x 10 m /s 
(15 SCFM) compressor, the prototype nozzle was designed for that capacity 


TABLE 4. Corona Scaling Laws. 


Increase all dimensions by a factor a: 

^e ,a a 

I. = al. 

J J 

Increase pressure from Pj^ to p: 

Eu = E. breakdown field (Paschen's law) 

b b,N Puj 


^e ■ ^e,N ^ 


V = V 

NPn 

I = I 


potential on attractor 
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with an exit Mach number of 1.35. The other expected characteristics of 
the prototype nozzle are given in Table 2. 

The unit has been designed to allow for the variation of the fea- 
tures listed in Table 1. In addition, attention has been paid to making 
the unit easily modified to allow optimization of the design. 
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3.0 DRAWINGS OF PROTOTYPE CHARGED PARTICLE FOG DISPERSAL UNIT 


Detailed drawings of the prototype unit are shown in Figures 9, 10, 
and 11. The nozzle unit will be constructed of machined and polished 
Plexiglas units that will be sealed with 0-rings. The unit will be 
bolted together using threaded rods. Connections to the prototype 
nozzle will be constructed from conmercially available pipe fittings. 

The entire control unit and the nozzle will be mounted on handtrucks to 
make the unit portable. 


20 



Figure 9 Cross section of prototype nozzle. 


NOTE: All dimensions in inches. 



Figure 10 Details of prototype nozzle. 
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Figure 11 Details of prototype fog dispersal system. 





4.0 DESIGN OF VERIFICATION EXPERIMENTS FOR FOG DISPERSAL UNIT 
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exit [19]. The particle size will be inferred from measured charge to 
mass ratios and by using diffusion charging or electric field breakdown 
theory. Also the possibility of droplet size measurement with laser 
techniques will be investigated. The attractor, needle, and nozzle exit 
currents will be measured with commercial microammeters. In addition. 



TABLE 5. Quantities to be Measured. 


Nozzle Properties: 

• Stagnation pressure 

t Stagnation temperature 

• Stagnation relative humidity 

• Exit Mach number 

• Exit droplet size 

Electrical Properties: 

• ^C’ ^A’ 

• Charge density distribution 

• Electric field distribution 


Atmospheric Properties: 

• Pressure 

• Temperature 

• Relative humidity 

• Visibility 

• Wind velocity 


the electric field and the charge density external to the nozzle must be 
measured before the nozzle is turned on and during the operation of the 
nozzle. Ideally, the distribution of both quantities must be measured 
to verify the various theories about the external field [2]. This is 
difficult to do because of the interference between the measuring instru- 
ment and its support and connection to the ground and the field to be 
measured, but these effects can be minimized by calibration of the 
instruments in place. 

The vertical component of the static electric field will be measured 
using field mills [32,33,34,35]. The number of field mills to be used 
will depend somewhat on their availability. The field mills will be 

C 

required to either have a large dynamic range (10 to 10 volt/m) or two 
different sets of instruments will have to be used so that both the 
background field and the field due to the charged droplets can be mea- 
sured. The field mills will have to be calibrated both in the laboratory 
and in place in the field. They will be calibrated by fitting large 
plates to them and applying a DC voltage between the plates. Comparison 
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with a laboratory-calibrated field mill will be required in the field 
because of the influence of the field mounting upon the calibration. 

It is also necessary to measure the space-charge distribution in 
the atmosphere in the region of the nozzle. A cage technique was used 
by Gourdine [5] and appears to be the most useful method to measure the 
charge density at considerable heights above the ground [36]. Calibra- 
tion is difficult but comparison can be made between the estimates of 
the charge density using the measurements of the electrical field and 
Poisson's equation and the cage measurement. 

The remaining atmospheric measurements indicated in Table 5 will be 
made using standard techniques. The complete field program was described 
previously [2]. 


26 


5.0 CONCLUDING REMARKS 


The detailed design of a prototype charged particle fog dispersal 
unit has been presented in this report. The concepts used for the 
design are well understood and have been widely used in the past both 
for fog dispersal units and for electrofluid dynamic power generation. 
Therefore, it is anticipated that the prototype unit will work effi- 
ciently as a current source over a wide range of conditions. However, 
the characteristics of the prototype unit that will ensure efficient 
dispersal of an atmospheric fog of given properties are still not known. 
In particular, the optimum water droplet size is unknown. Part of the 
uncertainty results from a lack of knowledge of how a fog is dispersed 
by the charged droplets. It is possible that it is dispersed by coales- 
cence and subsequent movement in the gravitational field to the ground 
[9,10] or the droplets could transfer part of their charge to the fog 
particles causing them to follow the electric field to the ground [30]. 
Also, it is possible that only some coalescence is required to increase 
the visibility through the fog without precipitation. Different mecha- 
nisms quite possibly could be involved in fog dispersal in different 
types of fogs. These questions must be addressed in a full field test 
of a group of prototype units. Only a sufficiently large group of units 
which create a nearly uniform electric field distribution in the hori- 
zontal direction will fully verify the two-dimensional fog dispersal 
theories [30]. The number of units to be used is dependent on the 
performance results from the prototype unit. Thus, the more complete 
field tests should be performed only after complete examination of the 
single unit that is described in Section 4.0. 
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